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The Green's functions of QCD encode important information about the infrared dynamics of the theory. The 
main non-perturbative tools used to study them are their own equations of motion, known as Schwinger-Dyson 
equations, and large-volume lattice simulations. We have now reached a point where the interplay between 
these two methods can be most fruitful. Indeed, the quality of the lattice data is steadily improving, while 
a recently introduced truncation scheme for the Schwinger-Dyson equations makes their predictions far more 
reliable. In this talk several of the above points will be reviewed, with particular emphasis on how to enforce the 
crucial requirement of gauge invariance at the level of the Schwinger-Dyson equations, the detailed mechanism of 
dynamical gluon mass generation and its implications for the ghost sector, the non-perturbative effective charge 
of QCD, and the indirect extraction of the Kugo-Ojima function from existing lattice data. 



1. Introduction 

The basic building blocks of QCD are the 
Green's (correlation) functions of the funda- 
mental physical degrees of freedom, gluons and 
quarks, and of the unphysical ghosts. Even 
though it is well-known that these quantities are 
not physical, since they depend on the gauge- 
fixing scheme and the parameters used to renor- 
malize them, it is widely believed that reliable 
information on their non-perturbative structure 
is essential for unraveling the infrared dynamics 
of QCD. Indeed, in addition to their relevance 
for phenomenology, the QCD Green's function 
encode information on confinement, albeit in a 
rather subtle way. 

The two basic non-perturbative tools that per- 
mit the exploration of the infrared domain of 
QCD are (i) the lattice, where space-time is dis- 
cretized and the quantities of interest are evalu- 
ated numerically, and (ii) the infinite set of in- 
tegral equations governing the dynamics of the 
QCD Green's functions, known as Schwinger- 
Dyson equations (SDE). Given that both the lat- 
tice and the SDE aspire to describe essentially 
the same physics, it is important to advance their 
complementarity and strengthen their mutual in- 



terplay. In fact, it would seem that we have 
reached a point in time where the meaningful 
and systematic comparison between lattice and 
SDE results constitutes a tangible reality. Indeed, 
the quality of lattice data is steadily improving, 
while, due to several recent developments p], we 
have at our disposal, for the first time, a mani- 
festly gauge invariant truncation scheme for the 
SD series of QCD. 

It it generally accepted by now that the lattice 
yields in the Landau gauge (LG) an infrared fi- 
nite gluon propagator and an infrared finite (non- 
enhanced) ghost dressing function. This rather 
characteristic behavior has been firmly estab- 
lished recently using large-volume lattices, for 
pure Yang-Mills (no quarks included), for both 
SU(2) [2] and SU(3) [3]. In this talk we will 
review the SD part of this story, within the 
truncation scheme based on the pinch technique 
(PT) [4I5I6| and its connection with the back- 
ground field method (BFM) [7]. We discuss the 
plethora of conceptual issues and the wide range 
of possibilities that emerge when the dynamical 
gluon mass generation picture of QCD is adopted. 
Moreover, we clarify some basic but subtle is- 
sues related with the non-perturbative definition 
of the QCD effective charge. 
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2. Problem(s) with the conventional SDEs 

The SDEs provide a formal framework for 
tackling physics problems requiring a non- 
perturbative treatment. Given that the SDEs 
constitute an infinite system of coupled non-linear 
integral equations for all Green's functions of the 
theory, their practical usefulness hinges crucially 
on one's ability to devise a self-consistent trun- 
cation scheme that would select a tractable sub- 
set of these equations, without compromising the 
physics one hopes to describe. Devising such a 
scheme, however, is very challenging, especially in 
the context of non-Abelian gauge theories, such 
as QCD. For the purposes of this presentation we 
will only focus on the aspect of the problem re- 
lated with gauge-invariance 

In Abelian gauge theories the Green's functions 
satisfy naive Ward identities (Wis) : the all-order 
generalization of a tree-level WI is obtained by 
simply replacing the Green's functions appearing 
in it by their all-order expressions. In general, 
this is not true in non-Abelian gauge theories, 
where the Wis arc modified non-trivially beyond 
tree-level, and are replaced by more complicated 
expressions known as Slavnov- Taylor identities 
(STIs): in addition to the original Green's func- 
tions appearing at tree-level, they involve various 
composite ghost operators. 

To appreciate why the Wis are instrumental for 
the consistent truncation of the SDEs, while the 
STIs complicate it, let us first consider how nicely 
things work in an Abelian case, namely scalar 
QED (photon), and then turn to the complica- 
tions encountered in QCD (gluon). 

The full photon (gluon) propagator, in a Re- 
type of gauge, is given by 



(i) 



where P^ v (q) — g^ v — q^q v /q 2 is the transverse 
projector. The scalar function A(q 2 ) is related to 
the all-order (photon) gluon self-energy n M „ (q) = 
P^{q)Ii{q 2 ) through A" 1 {q 2 ) = q 2 + lU{q 2 ). We 
also define the dimensionlcss vacuum-polarization 
U{q 2 ) 7 as IL(q 2 ) = q 2 U(q 2 ). 

Note now a crucial point: local gauge- 
invariance (BRST in the case of the gluon) forces 
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Figure 1. The SDE for the photon self-energy in 
scalar QED. 



n MV (q) (photon and gluon alike) to satisfy the 



fundamental transversality relation 
qm^{q) =0, 



(2) 



both perturbatively (to all orders) and non- 
perturbatively. 

The SDE governing H lll ^(q) in scalar QED is 
shown in Fig[T] The main question we want to ad- 
dress is the following: can one truncate the rhs of 
FiglU i.e. eliminate some of the graphs, without 
compromising the transversality of (q) ? The 
answer is shown already in FigQ] the two blocks 
of graphs [(ai) + (a 2 )] and [(a 3 ) + (o 4 ) + (a 5 )] are 
individually transverse, i.e. 
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0. 



(3) 



The reason for this special property are precisely 
the naive Wis satisfied by the full vertices ap- 
pearing on the rhs of the SDE; for example, the 
first block is transverse simply because the full 
photon-scalar vertex [white blob in (ai)] sat- 
isfies the WI 



q^T fl = e[V- 1 (k + q)-V- 1 (k)], 



(4) 



where V(q) is the full propagator of the charged 
scalar. A similar WI relating the four-vertex with 
a linear combination of forces the transversal- 
ity of the second block in Fig[TJ Thus, due to the 
simple Wis satisfied by the vertices appearing on 
the SDE for H^(q), one may omit the second 
block of graphs and still maintain the transver- 
sality of the answer intact, i.e. the approximate 
n M „(<7) obtained after this truncation satisfies 
Let us now turn to the conventional SDE for 
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Figure 2. The conventional SDE for the gluon 
self-energy 
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Figure 3. Diagrammatic representation of the 
functions A and H. 



the gluon self-energy, in the gauges, given in 
FiglU Clearly, by virtue of ([2]), we must have 



(5) 



i=l 



However, unlike the Abelian example, the dia- 
grammatic verification of i.e. through con- 
traction of the individual graphs by q^, is very 
difficult, essentially due to the complicated STIs 
satisfied by the vertices involved. For example, 
the full three-gluon vertex T allu (q, k\, k%) satisfies 
the STI 

q a r a ^ = F(q)A- 1 (kf)PJ(k 2 )H^(k u k 2 ) 

- - F{q)A- 1 {k 2 1 )P;{k 1 )H vl {k 2 , h), (6) 

where F(q) = q 2 D(q) is the ghost dressing func- 
tion, D(q) is the ghost propagator, and H^ u , de- 
fined in FigGH is related to the full gluon-ghost 
vertex by q v H^ v (k, q) — ~iT„(k,q); at tree-level, 
H^J = ig^v In addition, some of the perti- 
nent STIs are either too complicated, such as 
that of the conventional four-gluon vertex, or 
they cannot be cast in a particularly convenient 
form, such as the STI of the conventional gluon- 
ghost vertex. The main practical consequence 
of all this is that one cannot truncate the rhs 
of Fig(2] in any obvious way without violating 
the transversality of the resulting 11^ (g). For 
example, keeping only graphs (ai) and {a 2 ) is 
not correct even perturbatively, since the ghost 
loop is crucial for the transversality of n M „ al- 
ready at one-loop; adding (03) is still not suffi- 
cient for a SDE analysis, because (beyond one- 
loop) 9"[(ai) + (a 2 ) + (a 3 )W + 0. 



3. Truncating gauge-invariantly: 
The PT-BFM framework 

Recently a new truncation scheme for the SDEs 
of QCD has been proposed, that respects gauge 
invariance at every level of the "dressed-loop" ex- 
pansion. This becomes possible due to the drastic 
modifications implemented to the building blocks 
of the SD series, i.e. the off-shell Green's func- 
tions themselves, following the general methodol- 
ogy of the PT [4I5I6] , The PT is a well-defined 
algorithm that exploits systematically the BRST 
symmetry in order to construct new Green's func- 
tions endowed with very special properties. In 
particular, the crucial property for our discussion 
is that they satisfy Abelian Wis instead of the 
usual STIs The PT may be used to rearrange sys- 
tematically the entire SD series P. In the case 
of the gluon self-energy it gives rise dynamically 
to a new SDE, Fig|3J similar to that of FiglH but 
with two important differences: 

(i) What appears on the lhs is not the conven- 
tional self-energy 11^, but rather the PT-BFM 
self-energy, denoted by ILjj,. Of course, IT M „ is 
also transverse, i.e. q M II Ml/ (g) = 0. Quite inter- 
estingly, the two quantities are connected by a 
powerful formal identity [5] stating that 

A(q 2 ) = [l + G(<z 2 )] 2 A(g 2 ), (7) 
with G(q 2 ) defined from A /lu (q), FigOl 

iA„ v (q)=X f H^D(k + q)AP-(k)H (TV (k,q), 

Jk 



(8) 
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where A = g 2 C^, with Ca the Casimir eigen- 
value of the adjoint representation [Ca = N for 
SU(N)], and J k = (i 2s (2Tr)- d Jd d k, with d = 4-e 
the dimension of space-time. 

(ii) The graphs appearing on the rhs con- 
tain the conventional self-energy ILj„ as before, 
but are composed out of new vertices (Figd]). 
These new vertices correspond precisely to the 
Feynman rules of the BFM [7], i.e., it is as if 
the external gluon had been converted dynam- 
ically into a background gluon. As a result, 
the full vertices f^{q,kj,k 2 ), T a a nm (q, k u k 2 ), 

r^;(9,fci,fe,fc 3 ), and r— fa,**,*!,, fe 3 ) ap- 
pearing on the rhs of the SDE shown in Fig[4] 
satisfy the simple Wis 

9 a fr = i5/ aTOn [i>- 1 (Ai)-i3" 1 (fe)]. 

<fK7Z = gf adr ^{d + k2,k 3 M) + ^, 

qaf arnnr = fl /«^r^( ? + ft, A*, + cp. (9) 

where "cp" stands for "cyclic permutations" The 
implication of this for the truncation of the SDE 
are far-reaching. Indeed, using these Wis, it is 
elementary to show that pQ 



{ai) flu = ) q^y^^ai)^ = . 
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(10) 



i=7 



Evidently, the SDE is composed of "one-loop" 
and "two-loop" dressed blocks that are indi- 
vidually transverse, exactly as happened in the 
Abelian case. In fact, the resulting pattern is even 
more spectacular: the gluon and ghost diagrams 
form separate transverse blocks! 

4. Gluon mass generation: the big picture 

Even though the gluon is massless at the level 
of the fundamental QCD Lagrangian, and re- 
mains massless to all order in perturbation the- 
ory, the non-perturbative QCD dynamics gen- 
erate an effective, momentum-dependent mass, 
without affecting the local SU(3) C invariance, 
which remains intact [1]. 
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Figure 4. The SDE for the gluon self-energy in 
the PT-BFM framework. 



The gluon mass generation is a purely non- 
perturbative effect. At the level of the SDEs the 
generation of such a mass is associated with the 
existence of infrared finite solutions for the gluon 
propagator |4l9j . i.e. solutions with A~ 1 (0) > . 
Such solutions may be fitted by "massive" propa- 
gators of the form A -1 (q 2 ) = q 2 +m 2 (q 2 ): m 2 (q 2 ) 
is not "hard", but depends non-trivially on the 
momentum transfer q 2 . In addition, one obtains 
the non-perturbative generalization of a(q 2 ), the 
QCD running coupling (effective charge), in the 
form 4 



a- 1 {q 2 )=b\n 



4m 2 (q 2 



A 2 



(11) 



The m 2 (q 2 ) in the argument of the log- 
arithm tames the Landau pole, and a{q 2 ) 
freezes at a finite value in the IR, namely 
O! -1 (0) = 61n(4m 2 (0)/A 2 ). Moreover, the gluon 
mass forces F(q 2 ) to stay IR-finite (non- 
enhanced) [9TTU] . 

(i) The Schwinger mechanism 

In order to obtain massive solutions gauge- 
invariantly, it is necessary to invoke the well- 
known Schwinger mechanism |llj : if, for some 
reason, II (q 2 ) acquires a pole at zero momen- 
tum transfer, then the vector meson becomes 
massive, even if the gauge symmetry forbids a 
mass at the level of the fundamental Lagrangian. 
Indeed, it is clear that if the vacuum polariza- 
tion U(q 2 ) has a pole at q 2 = with positive 
residue /i 2 , i.e. Tl(q 2 ) = fi 2 /q 2 , then (in Eu- 
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Figure 5. Vertex with non-perturbative massless 
excitations triggering the Schwinger mechanism. 



clidean space) A~ 1 (<7 2 ) = q 2 + [i 2 . Thus, the 
vector meson becomes massive, A _1 (0) = [i 2 , 
even though it is massless in the absence of in- 
teractions (g = 0). There is no physical principle 
which would preclude Tl(q 2 ) from acquiring such 
a pole. In a strongly- coupled theory like QCD this 
may happen for purely dynamical reasons, since 
strong binding may generate zero-mass bound- 
state excitations [12] . The latter act like dynam- 
ical Nambu-Goldstone bosons, in the sense that 
they are massless, composite, and longitudinally 
coupled; but, at the same time, they differ from 
Nambu-Goldstone bosons as far as their origin is 
concerned: they do not originate from the spon- 
taneous breaking of any global symmetry [3] . 

As we will see in a moment, the exact way 
how the Schwinger mechanism is integrated into 
the SDE is through the form of the three-gluon 
vertex. In particular, one assumes that the ver- 
tex contains dynamical poles ~ l/q 2 [see FigJS], 
which will trigger the Schwinger mechanism when 
inserted into the SDE for the gluon self-energy 
(other vertices, e.g, the four-gluon vertex, may 
have similar poles). 

(ii) Gluon mass and confinement 

It has been occasionally argued that the con- 
cept of a massive gluon may be at odds with con- 
finement, because a massive gauge field does not 
give rise to the long-range potential necessary for 
this latter phenomenon to occur. This is, how- 
ever, not true: the gluon mass does not obstruct 
confinement, it enables it! 

Of course, the exact way how this happens is 
very intricate, and is inextricably connected with 
the notion of a quantum soliton. A quantum soli- 
ton is a localized finite-energy configuration of 
gauge potentials arising from an effective action 



- c * } ; = (- — ► — r x + --►-^-0"»-^-- 

q q q k+q 

Figure 6. The SDE for the ghost propagator. 

that summarizes quantum effects not present in 
the classical action; in our case, the quantum ef- 
fect is a gauge-invariant dynamical gluon mass. 
Specifically, an effective low-energy field theory 
for describing the gluon mass is the gauged non- 
linear sigma model known as "massive gauge- 
invariant Yang-Mills" [T3j , with Lagrangian den- 
sity C m 

Cm = \Gl v -m 2 Ti [A, - g- 1 U(9)d l± U- 1 (9)] 2 , (12) 

where = ^ £ a X a A a ^ the A a are the SU(3) 
generators (with TrA a Ab = 26 a b), and the N x 
N unitary matrix U(9) = exp [iiA a 6> a ] describes 
the scalar fields 9 a . C m is locally gauge-invariant 
under the combined gauge transformation 

U' = U{9') = VU{9), (13) 

for any group matrix V = exp [iiA a w a (x)] , 
where Lo a (x) are the group parameters. C m ad- 
mits vortex solutions - these are the aforemen- 
tioned quantum solitons - with a long-range pure 
gauge term in their potentials, which endows 
them with a topological quantum number corre- 
sponding to the center of the gauge group [Zn 
for SU(N)], and is, in turn, responsible for quark 
confinement and gluon screening. Specifically, 
center vortices of thickness ~ m _1 , where m is 
the induced mass of the gluon, form a conden- 
sate because their entropy (per unit size) is larger 
than their action. This condensation furnishes an 
area law to the fundamental representation Wil- 
son loop, thus confining quarks [13] . 

5. Solutions of the SDEs and comparison 
with the lattice results 

We are now in a position to incorporate into 
the SDE truncation scheme based on the PT- 
BFM, described above, the concrete dynamics 
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Figure 7. Comparison between SDE results of [9] and the lattice results of [3]. 



that will give rise to an infrared finite gluon prop- 
agator, signaling the dynamical generation of a 
gluon mass. To that end, we truncate (gauge- 
invariantly!) the SDE for the gluon propagator 
by keeping only the one-loop dressed contribu- 
tions, i.e. the first two blocks of graphs in FigJU 
Then, we have [9] 



a-V)Pm«(?) = 



[1 + G(<f)]2 



(14) 



and 



We next express r Ma| g(g, k\, k 2 ) 
T /1 (q,ki,k 2 ) [appearing in (oi) and (a 3 ) of FigJH 
respectively] as a function of the gluon and ghost 
self-energy, respectively, in such a way as to 
automatically satisfy the first two Wis of J9}; 
failure to do so would invariably compromise the 
transversality of the answer. The Ansatz we will 
use for T ain ,(q, ki,k 2 ) is 



r^ = r^ + i^[n^(fc 2 ) 



iMM] - (is) 



and a similar expression for r„. The essential 
feature is the presence of massless pole terms, re- 
quired for triggering the Schwingcr mechanism. 



The resulting final expression for the SDE of (Tj"4"|) 
is too lengthy to report here [9]. The main point 
is that, as desired, A _1 (0) > 0, i.e. the gluon 
self-energy is infrared finite. 

For the conventional ghost-gluon vertex T^, ap- 
pearing in the ghost SDE of Fig[51 we use its tree- 
level expression, i.e., = — p M ; this is perfectly 
legitimate, since injhe PT-BFM formalism the 
two ghost vertices, and r„, are different. Fi- 
nally, for we use its tree-level value, ig^ v . 

Then, we obtain (Euclidean space) [9115] 



F-V) 
G(q 2 ) 
Li/) 



1 - A 

3 A 
3 A 



(g-fc) 2 

q 2 k 2 

(q ■ kf 
k 2 q 2 

Aq-k)'- 



k 2 q 2 



A(k)D(q + k), 

A(k)D(k + q), 

A(k)D(k + q). 

(16) 



In fact, there exists a powerful formal identity 
relating F(q 2 ), G(q 2 ), and L(q 2 ), namely [14] 



F- 1 (q 2 ) = l + G(q 2 ) + L( q 2 ). 



(17) 
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In addition to its formal derivation [T3] , the above 
relation has been recently obtained at the level 
of the SDEs defining these three quantities [TS] . 
Adding by parts Eqs. (|16p we can verify that 
Eq. (|17[) is indeed satisfied under the approxima- 
tions employed. 

Eq. (|17| merits further analysis. The origin of 
this equation is the BRST symmetry of the the- 
ory [T3]; in that sense, Eq. ([T7|) has the same na- 
ture as the STIs. Therefore, just as happens with 
the latter, (fl7|) must not get deformed after renor- 
malization [15]. Thus, denoting by Z u the renor- 
malization constant relating the bare and renor- 
malized functions, Aq 1 ' and A M1/ , [see ((8])] through 

Z u [g^ + A^(q)]=g^+A^(q), (18) 

the requirement of non-deformation im- 
poses the crucial condition Z u = Z c [IB"] . 

Thus, Z c [l + G (q 2 ,A 2 )} = 1 + G(? 2 ,M 2 ) and 
Z c L (q 2 ,A 2 )^L(q 2 ^ 2 ). 

The solutions obtained form the above system 
of SDEs is shown in Fig[71 where it is compared 
with the corresponding lattice data [L(q 2 ) is nu- 
merically suppressed and is not shown here] . Note 
that while there is good qualitative agreement 
with the lattice, there is a significant discrepancy 
(a factor of 2) in the intermediate region of mo- 
menta. This of course may not come as a sur- 
prise, given that the "two-loop" dressed part of 
the SDE for A has been omitted [the last two 
blocks in FigfJ]. Even though this omission has 
not introduced artifacts (since it was done gauge- 
invariantly), the terms left out are expected to 
modify precisely the intermediate region, given 
that both the IR and UV limits of the solu- 
tions are already captured by the terms consid- 
ered here. 

6. The effective charge of QCD 

Due to the Abelian Wis satisfied by the PT- 
BFM Green's functions, the new A~ 1 (q 2 ) absorbs 
all the renormalization-group (RG) logarithms, 
exactly as happens in QED with the photon self- 
energy [4 5 6. As a result, the renormalization 
constants of the gauge-coupling and of the PT 
gluon sclf-cncrgy, defined as 

gill 2 ) = Z- 1 ([i 2 )g , 



A(q 2 ^ 2 ) = ^V)A (g 2 ), 



(19) 



where the "0" subscript indicates bare quanti- 

^ — 1/2 

ties, satisfy the QED-likc relation Z g = Z A ' . 
Thus, regardless of the renormalization prescrip- 
tion chosen, the product 

d (Q 2 ) = dlMl 2 ) = 3 V)A(q 2 , M 2 ) = d(q 2 )(20) 

retains the same form before and after renor- 
malization, i.e., it forms a RG-invariant (fi- 
independent) quantity [4j. Note that, by virtue 
of ([7]), we have that 



d(q 2 )=g 2 ^ 2 ) 



A(g 2 ,/x 2 



[l + G(g 2 , M 2 



(21) 



To see how beautifully the /i-independence of 
d(q 2 ) is captured by the solutions of the SDEs 
we study, note in Fig[5]the explicit /i-dependence 
of the individual ingredients entering into its def- 
inition. However, when they are put together ac- 
cording to the rhs of (f2"Tj) . the resulting d(q 2 ) is 
practically /i-independent! 

The next step is to extract out of d(q 2 ) a di- 
mensionless quantity, that would correspond to 
the QCD effective charge. Pcrturbatively, i.e. for 
asymptotically large momenta, it is clear that the 
mass scale is saturated simply by q 2 , the bare 
gluon propagator, and the effective charge is de- 
fined by pulling a 1/q 2 out of the corresponding 
RG-invariant quantity, according to 



d(q 2 ) 



til 2 



(22) 



where g 2 (q 2 ) is the RG-invariant effective charge 
of QCD; at one-loop 



sV) 



1 



l + bgnn(qy^) 61n(g 2 /A 2 CD ) 



.(23) 



where A QCD denotes an RG-invariant mass scale 
of a few hundred MeV. 

Of course, given that in the IR the gluon propa- 
gator becomes effectively massive, particular care 
is needed in deciding exactly what combination of 
mass-scales ought to be pulled out. It would cer- 
tainly be unwise, for example, to continue defin- 
ing the effective charge by forcing out just a fac- 
tor of l/<7 2 ; such a procedure would furnish (triv- 
ially) a completely unphysical coupling, vanishing 
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in the deep infrared, where QCD is supposed to 
be strongly coupled ! The correct procedure, in- 
stead, is to factor out a "massive" propagator, of 
the form [q 2 + m 2 (q 2 )]^ 1 , i.e. one must set [4115] 

d V) = 2 tK] 2, > ( 24 ) 

where m 2 (q 2 ) is a momentum-dependent gluon 
mass. Clearly, for q 2 3> m 2 (q 2 ) the expression on 
the rhs of fl24|) goes over to that of (|22j). The 
effective charge, a PT (<7 2 ) = g 2 (g 2 )/47r, obtained 
from the d(g 2 ) after using is also shown in 
FigJSJ We have assumed "power-law" running of 



the mass, 

mV) = i^7> (25) 
q z + rag 

consistent with various independent studies ,16]. 

Let us finally comment on some important con- 
ceptual issues. As already mentioned, the "freez- 
ing" of the coupling is a direct consequence of the 
appearance of a mass in the RG logarithm. This 
fundamental property of the strong coupling may 
be reformulated in terms of what in the language 
of the effective field theories is referred to as "de- 
coupling" : at energies sufficiently inferior to their 
masses, the particles appearing in the loops (in 
this case the gauge bosons) cease to contribute to 
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the "running" of the coupling. Note, however, a 
crucial point, which has been the source of con- 
siderable confusion in the recent literature: the 
"decoupling" , as described above, is not a syn- 
onym for non-interactive! In the electroweak sec- 
tor, for example, such a "decoupling" takes place 
indeed, since below the mass of the W the gauge 
boson loops do not contribute to the running. 
However, this is by no means the same as saying 
that the theory is non-interactive (for one thing, 
the (3 decay still takes place.) Another central 
point is that when the QCD charge is constant 
(non-vanishing!) in the infrared (and the quark 
masses are ignored), QCD becomes conformally 
invariant, and the AdS/CFT correspondence be- 
comes applicable [17] . In conclusion: The gluon 
mass keeps QCD strongly coupled, and with the 
"conformal window" open! 

6.1. A case of SDE-lattice synergy: 
The Kugo-Ojima function. 

The Kugo-Ojima (KO) scenario [18] claims to 
establish a highly non-trivial link between con- 
finement and the infrared behavior of the ghost 
dressing function F{q 2 ). In a nutshell, a suffi- 
cient condition for the realization of the KO con- 
finement ( "quartet" ) mechanism is that a certain 
correlation function, u(q ), defined as 

J d A x e- l ^-^(T[(2V)™ {f nrs A n v c s ) y )) = 

P^{q)5 mn u{q 2 ), (26) 

should satisfy the condition u(0) = —1. Given 
that u(0) is related to the infrared behavior of 
F(q 2 ) through the identity F _1 (0) = 1 + u(0), 
the KO confinement scenario predicts a divergent 
ghost dressing function, and vice-versa. As was 
already mentioned, however, the ghost dressing 
function is not enhanced on the lattice [2|3 . In 
addition, and perhaps not surprisingly, the lat- 
tice [TH] finds no evidence of u(0) = —1 either: 
u{q 2 ) saturates in the deep infrared around ap- 
proximately —0.6. 

Quite remarkably, even though their field- 
theoretic origin appears to be completely differ- 
ent, the KO function coincides (in the LG!) with 
the G(q 2 ), introduced in ©, namely [14] 

u(q 2 ) = G{q 2 ). (27) 



One may then substitute into the SDE for G(q 2 ), 
given in f|16[>. the available lattice data on the 
gluon and ghost propagators (e.g. [3]), and ob- 
tain an indirect determination for u(q 2 ) [20j . Note 
that from ([T6]) one can verify explicitly [15] that 
L(0) = 0. The u(q 2 ) so obtained may be then 
compared with the direct lattice determination of 
the same quantity [19] . from its defining Eq. ([26]) . 

The outcome of this procedure is displayed in 
the left panel of FigOll Evidently, the coincidence 
between the result obtained from the combined 
(indirect) SDE-lattice appraoch and that of the 
pure (direct) lattice is rather good. Note in ad- 
dition an important (and only recently appreci- 
ated) point: u{q 2 ), and in particular u(0), de- 
pend explicitly on the renormalization point /it, as 
one would expect, given that the KO function is 
not an intrinsically /^-independent quantity [un- 
like, e.g., d(q 2 )]. The implications of this, and 
other facts, for the KO confinement mechanism 
(and approaches relying on it) are currently un- 
der intense scrutiny [2 1] . 

7. Conclusions 

In this presentation we have outlined the salient 
features of the SDEs formulated within the PT- 
BFM framework, and have given several examples 
of the considerable potential offered by their in- 
terplay with the lattice simulations. Clearly sev- 
eral aspects need to be further investigated, e.g. 

(i) the dependence of the infrared behavior on the 
gauge chosen (e.g., LG vs Feynman gauge); 

(ii) study on the lattice the auxiliary (ghost) 
Green's functions appearing in the new SDE; 
(Hi) refinement of the SDE treatment to improve 
the agreement with lattice; 

(iv) the possible connection between the SDEs 
and the Gribov horizon. 

We hope that the ongoing effort of the lattice and 
SDE communities will soon shed light on these 
and many more questions. 
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